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Abstract 

Although standard economics textbooks are seldom interested in 
production networks, modern economies are more and more based 
upon suppliers/customers interactions. One can consider entire sec- 
tors of the economy as generalised supply chains. We will take this 
view in the present paper and study under which conditions local 
failures to produce or simply to deliver can result in avalanches of 
shortage and bankruptcies across the network. We will show that a 
large class of models exhibit scale free distributions of production and 
wealth among firms and that metastable regions of high production 
are highly localised. 

1 Networks of firms 

Firms are not simply independent agents competing for customers 
on markets. Their activity involves many interactions, and some of 
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them even involve some kind of cooperation. Interactions among firms 
might include: 

• information exchange [HI HI, Eli 

• loans[7j,|Sl; 

• common endeavours 9 ; 

• partial ownership 0; 

• and of course economic transactions allowing production IQj (the 
present paper). 

Economic activity can be seen as occurring on an economic network 
("le tissu economique" ) : firms are represented by vertices and their 
interactions by edges. The edges are most often asymmetric (think for 
instance of providers/customers interactions). The availability of em- 
pirical data has provoked research on the structure of these networks: 
many papers discuss their "small world properties" and frequently 
report scale free distribution [2j of the connections among firms. 

The long term interest of economic network research is rather the 
dynamics creating or occurring on these nets: how are connections 
evolving, what are the fluxes of information, decisions economic 
transactions etc ... But dynamic studies lag behind statistical ap- 
proaches because of conceptual difficulties and because time series 
of individual transactions are harder to obtain than time aggregated 
statistics. 

The recent cascade of bankruptcies that occurred in Eastern Asia 
in 1997, provoked some research on the influence of the loans network 
structure on the propagation of "bad debts" and resulting avalanches 
of bankruptcies ([Zj,[Hj) • One of the most early papers on avalanche 
distribution in economic networks is due to Bak et al |1U| . It con- 
cerns production networks: edges represent suppliers/customers con- 
nections among firms engaged in batch production activity. The au- 
thors describe the distribution of production avalanches triggered by 
random independent demand events at the output boundary of the 
production network. 

These papers ([ZMBl and[TU]) are not based on any empirical de- 
scription of the network structure, but assume a very simple interac- 
tion structure: star structure in the first casejZj,[H], periodic lattice in 
Bak et al paper ^U]. They neither take into account price dynamics. 

The present paper is along these lines: we start from a very sim- 
ple lattice structure and we study the consequences of simple lo- 
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cal processes of orders/production (with or without failure) /delivery/ 
profit /investment on the global dynamics: evolution of global produc- 
tion and wealth in connection to their distribution and local patterns. 
In the spirit of complex systems analysis, our aim is not to present spe- 
cific economic prediction, but primarily to concentrate on the generic 
properties (dynamical regimes, transitions, scaling laws) common to 
a large class of models of production networks. 

A minimal model of a production network will first be introduced 
in section 2. Simulation results are presented in section 3. Section 4 
is a discussion of the genericity of the obtained results: reference is 
made to comparable soluble models. We also summarise the results of 
several variants of the simplest model. The conclusion is a discussion 
of possible applications to geographical economics. 

2 A simple model of a production net- 
work 

We can schematise the suppliers/customers interactions among firms 
by a production network, where firms are located at the vertices and 
directed edges represent the delivery of products from one firm to its 
customers (see figure 1). 

Independent local failures to produce (or to deliver) by a firm might 
give rise to the propagation of shortage across the production network. 

We have chosen a simple periodic lattice with three input connec- 
tions of equal importance and three output per firm. The network is 
oriented from an input layer (say natural resources) towards an output 
layer (say the shelves of supermarkets). The transverse axis can be 
thought as representing either geographical position or some product 
space while the longitudinal axis relates to production. We here use 
a one dimensional transverse space to facilitate the representation of 
the dynamics by two-dimensional patterns, but there is no reason to 
suppose geographical or product space to be one-dimensional in the 
real world. 

In real economies, the network structure is more heterogenous with 
firms of unequal importance and connectivity. Furthermore some de- 
livery connections go backwards. Most often these backward connec- 
tions concern equipment goods; neglecting them as we do here implies 
considering equipment goods dynamics as much slower than consump- 
tion goods dynamics. Anyway, since these backward connections enter 
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Figure 1: Firms are located at the nodes of the lattice. Production (Y D ) 
flows from the resource input layer (A; = I) to the output layer (A; = 0), 
orders (Y) flow backward. 
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positive feedback loops, we have no reason to suppose that they would 
qualitatively disrupt the dynamics that we further describe. 

At each time step two opposite flows get across the lattice: orders 
are first transmitted upstream from the output layer; production is 
then transmitted downstream from the input layer to the output layer. 

• Orders at the output layer 

We suppose that orders are only limited by the production ca- 
pacity 1 Aoi of the firm in position 0, i, where indicates the 
output layer, and i the transverse position in the layer. 

Y 0i = q-Aoi (1) 

Yoi is the order in production units, and q a technological pro- 
portionality coefficient relating the quantity of product Y to the 
production capacity A, combining the effect of capital and labor. 
q is further taken equal to f without loss of generality. 

• Orders 

Firms at each layer k, including the output layer, transfer orders 
upstream to get products from layer k + 1 allowing them to pro- 
duce. These orders are evenly distributed across their 3 suppliers 
upstream. But any firm can only produce according to its own 
production capacity A^. The planned production is then a 
minimum between production capacity and orders coming from 
downstream: 

Y ki = min(q ■ Au,^ ^ k V>1 ) (2) 

v 6 

v stands for the supplied neighborhood, here supposed to be the 

three firms served by firm k, i (see figure 1). 

We suppose that resources at the input layer are always in excess 

and here too, production is limited only by orders and production 

capacity. 

• Production downstream 

Starting from the input layer, each firm then starts producing 
according to inputs and to its production capacity; but produc- 
tion itself is random, depending upon alea. We suppose that at 

1 A number of simplifying assumptions of our model are inspired from [5], especially 
the assumption that production is limited by production capacity, not by market. 
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each time step some catastrophic event might occur with con- 
stant probability V and completely destroy production. Fail- 
ures result in canceling production at the firm where they occur, 
but also reduce production downstream, since firms downstream 
have to reduce their own production by lack of input. These 
failures to produce are uncorrelated in time and location on the 
grid. Delivered production Y k \ by firm k, i then depends upon 
the production delivered upstream from its delivering neighbor- 
hood v\ at level k + 1: 



y& = (E y (W v Yki v >- € ^ (3) 

i>evl 2^i"e Vi , 1 ki" 

— Whenever any of the firms i' £ v[ at level k + 1 is not able to 
deliver according to the order it received, it delivers down- 
stream at level k to its delivery neighbourhood in pro- 
portion of the initial orders it received, which corresponds 
to the fraction term; 

— e(i) is a random term equals to with probability V and 1 
with probability 1 — V. 

The propagation of production deficit due to local independent 
catastrophic event is the collective phenomenon we are interested 
in. 

Profits and production capacity increase 

Production delivery results into payments without failure. For 
each firm, profits are the difference between the valued quantity 
of delivered products and production costs, minus capital decay. 
Profits Hfcj are then written: 

n kl = p-Y'-c-Y'-XAk (4) 

where p is the unit sale price, c is the unit cost of production, and 
A is the capital decay constant due to interest rates and material 
degradation. We suppose that all profits are re-invested into 
production. Production capacities of all firms are thus upgraded 
(or downgraded in case of negative profits) according to: 

A kl (t+1) = A kl (t)+U kl (t) (5) 
Bankruptcy and re-birth. 
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We suppose that firms which capital becomes negative go into 
bankruptcy. Their production capacity goes to zero and they 
neither produce nor deliver. In fact we even destroy firms which 
capital is under a minimum fraction of the average firm (typ- 
ically 1/500). A re-birth process occurs for the corresponding 
vertex after a latency period: re-birth is due to the creation of 
new firms which use the business opportunity to produce for the 
downstream neighborhood of the previously bankrupted firm. 
New firms are created at a unique capital, a small fraction of the 
average firm capital (typically 1/250). 2 . 

The dynamical system that we defined here belongs to a large class 
of non linear systems called reaction-diffusion systems (see e.g. |18j ) 
from chemical physics. The reaction part here is the autocatalytic loop 
of production and capital growth coupled with capital decay and death 
processes. The diffusion part is the diffusion of orders and production 
across the lattice. We can a priori expect a dynamical behaviour 
with spatio-temporal patterns, well characterised dynamical regimes 
separated in the parameter space by transitions or crossovers, and 
scale free distributions since the dynamics is essentially multiplicative 
and noisy. These expectations guided our choices of quantities to 
monitor during simulations. 

3 Simulation results 

3.1 Methods and parameter choice 

Unless otherwise stated, the following results were obtained for a pro- 
duction network with 1200 nodes and ten layers between the input 
and the output. 

Initial wealth is uniformly and randomly distributed among firms: 



One time step correspond to the double sweep of orders and pro- 
duction across the network, plus updating capital according to profits. 
The simulations were run for typically 5000 time steps. 

2 Adjusting these capital values relative to the average firm capital < A > is a standard 
hypothesis in many economic growth models: one supposes that in evolving economics 
such processes depend upon the actual state of the economv|lfi| and not upon fixed and 
predefined values. 




(6) 
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The figures further displayed correspond to: 

• a capital threshold for bankruptcy of < A > /500; 

• an initial capital level of new firms of < A > /250; 

Production costs c were 0.8 and capital decay rate A = 0.2. In 
the absence of failures, stability of the economy would be ensured 
by sales prices p = 1.0. In fact, only the relative difference between 
these parameters influences stability. But their relative magnitude 
with respect to the inverse delay between bankruptcy and creation of 
new firm also qualitatively influence the dynamics. 

In the limits of low probability of failures, when bankruptcies are 
absent, the linear relation between failure probability V and equilib- 
rium price p is written: 

p = c + \+ l --V (7) 

where I is the total number of layers. The | comes from the fact that 
the integrated damage due to an isolated failure is proportional to the 
average number of downstream layers. The slopes at the origin of the 
breakeven lines of figure 2 verify this equation. 

Most simulations were monitored online: we directly observed the 
evolution of the local patterns of wealth and production which our 
choice of a lattice topology made possible. Most of our understanding 
comes from these direct observations. But we can only display global 
dynamics or static patterns in this manuscript. 

3.2 Monitoring global economic performance 

The performance of the economic system under failures can be tested 
by checking which prices correspond to breakeven: the capital dynam- 
ics being essentially exponential, the parameter space is divided in two 
regions, where economic growth or collapse are observed. Drawing the 
breakeven manifolds for instance in the failure probability V and sale 
price p plane allows to compare the influence of other parameters . 
The growth regime is observed in the low V and high p region, the 
collapse regime in the high V and low p region. 

Figure 2 displays four breakeven manifolds corresponding to dif- 
ferent lattice depths. 

At low failure probability, the breakeven lines follow equation 7. 
At higher values of V, interactions among firms failures are important, 
hence the non linear increase of compensating prices. 
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Figure 2: Regime diagram in the sale price versus probability of failure plane. 
The time lag between bakruptcy and re-birth is 20. Two regions of growth 
and economical collapse at large times are separated by lines which position 
are fixed by simulation parameters. We here varied the production network 
depth: The red '+' line was obtained for a 3 layers net, the green 'x' line for 
a 5 layers net, the blue '*' line for an 8 layers net, and pink square line for a 
10 layers net. 
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Breakeven manifold are a simple test of the economic performances 
of the network: when performances are poor, the compensating sales 
price has to be larger. We checked for instance that increasing the 
bankruptcy threshold and new firms initial capital increase global eco- 
nomic performance. On the other hand, increasing the time lag be- 
tween bankruptcy and the apparition of new firms increase breakeven 
sale prices in the non-linear region. 

Among other systematic tests, we checked parent models with 
more realistic representations of production costs such as: 

• Influence of capital inertia; production costs don't instantly read- 
just to orders: capital and labour have some inertia which we 
modeled by writing that productions costs are a maximum func- 
tion of actual costs and costs at the previous period. 

• Influence of the cost of credit: production failures increase credit 
rates. 

Both variants of course yield higher breakeven sale prices; nevertheless 
these variants display the same generic properties that we will discuss 
in the next sections. 

Most further results, dynamical and statistical, are based on runs 
close to the breakeven price in order to avoid systematic drifts and 
recalibrations. 

3.3 Time evolution 

The simplest way to monitor the evolution of the system is to display 
the time variations of some of its global performance. Figure 3 displays 
the time variations of total delivered production Y d , total wealth A, 
total undelivered production due to failures and the fraction of active 
firms for a 1200x10 lattice, with a probability of failure of 0.05 and a 
compensation sale price of 1.185. Time lag between bankruptcy and 
and new firm creation is either 1 (for the left diagram) or 5 (for the 
right digram). 

The features that we here report are generic to most simulation at 
breakeven prices. During the initial steps of the simulation, here say 
1000, the wealth distribution widens due to the influences of failures. 
Bankruptcies cascades do not occur as observed by checking the num- 
ber of active firms, until the lowest wealth values reach the bankruptcy 
threshold. All quantities have smooth variations. Later, for t > 1000 
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Figure 3: Time evolution of wealth (red '+'), production (blue '*'), destroyed 
production (green 'x'), active firms (magenta empty squares) and production 
by the largest firm (cyan hollow squares). The network has 10 layers, 200 
firms per layer, V = 0.05 (the failure probability). The left diagram corre- 
sponds to a small time lag (1) between bankruptcy and firm re-birth, right 
diagram corresponds to a larger time lag (5). Vertical scale is logarithmic, 
which permits to have the four quantities displayed on the same time plot 
but reduces the apparent amplitude of fluctuations occurring when time is 
larger than 1000. 
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one observes large production and wealth fluctuations characteristic 
of critical systems. 

At larger time lag (5) between bankruptcy and firm re-birth, when 
bankruptcies become frequent, they can cascade across the lattice and 
propagate in both network directions as seen on the right diagram of 
figure 3. A surprising feature of the dynamics is that avalanches of 
bankruptcies are not correlated with production level. Even when only 
one tenth of the firms are active, the total production is still high. In 
fact, in this model, most of the total production is dominated by large 
firms, and avalanches which concern mostly small firms are of little 
consequence for the global economy. 

Battiston etal study more thoroughly the time dynamics of a re- 
lated model (large sale price fluctuations possibly inducing bankrupt- 
cies and lack of payment) in . 

3.4 Wealth and production patterns 

Like most reaction-diffusion systems, the dynamics is not uniform 
in space and display patterns. The wealth and production patterns 
displayed after 5000 time steps on figure 4 and 5 were obtained for 
V = 0.05 . They reflect wide distributions and spatial organisation. 
In these diagrams, production flows upward. The upper diagram dis- 
plays wealth A and the lower one production The intermediate 
bar is the colour scale, black=0, violet is the maximum wealth or pro- 
duction. (We in fact displayed square roots of A and Yd in order to 
increase the visual dynamics of the displays; otherwise large regions 
of the patterns would have been red because of the scale free distri- 
butions of A and Yd, see further). 

The important result is that although production has random fluc- 
tuations and diffuses across the lattice, the inherent multiplicative (or 
autocatalytic) process of production + re-investment coupled with lo- 
cal diffusion results in a strong metastable local organisation: the dy- 
namics clusters rich and productive firms in "active regions" separated 
by "poor regions" (in red or black). 

These patterns are evolving in time, but are metastable on a long 
time scale, say of the order of several 100 time steps as seen on the 
succession of production patterns at different steps of the simulation 
as one can observe on figure 6: successive patterns at time 1250, 1750 
and 2250. 

The relative importance of active (and richer) regions can be checked 
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Figure 4: Patterns of wealth(upper pattern) and production (lower pattern) 
after 5000 iterations steps with the parameter set-up of figure 3 (left) (time 
lag =1), for a 200x10 lattice. . For both patterns the output layer is the last 
one above. The intermediate line is the colour code, with minimal amplitude 
at the extreme left. We observe alternance of highly productive regions 
(in pink, blue and green colour), with less active regions (in red). Local 
production failures represented by black dots are randomly distributed across 
the production pattern. Only one bankrupted firm is observed on the wealth 
pattern. 
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Figure 5: Patterns of wealth(upper pattern) and production (lower pattern) 
after 5000 iterations steps with the parameter set-up of figure 3 (right) (time 
lag is 5). The same alternance of active and less active regions is observed, 
but with a larger time lag (5), we also get large zones of bankrupted firms in 
black. 




Figure 6: Successive patterns of wealth after 250, 750, 1250, 1750 and 2250 
time steps with the parameter set-up of figure 3 (right, time lag = 5) for a 
1200x10 lattice. 
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by a Zipf plot 12 . We first isolate active regions by "clipping" the dow- 
stream (along k axis) integrated wealth at a level of one thousandth 
of the total production 3 . 
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Figure 7: Separating richer regions. Downstream integrated wealth levels 
(green '+') are plotted as a function of their transverse position. The clipping 
level indicated by the red line isolates richer regions (those wealth peaks 
above the red line). 



We then transversally (along i axis) integrate the wealth of active 
regions and order these regional wealths to get the Zipf plots. 

All 3 Zipf plots display some resemblance with standard Zipf 12 
plots of individual wealth, firm size and city size. For the model 
discussed here, the size decrease following approximately a power law. 
The apparent 4 exponent is one when the time lag is 1. It is much 
higher when the time lag is 5. 

3 Clipping here means that when the production level is lower than the threshold it is 
set to zero 

4 the approximate algorithm that we use to isolate high productivity regions is respon- 
sible for the kinks in the Zipf plot 
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Figure 8: Zipf plot of wealth of the most active regions for the standard 
and adaptive firms models (cf. section 4.2). The vertical axis display the 
production relative to the total production. The red '+' correspond to the 
standard model with time lag = 5, green V to time lag = 1, and blue '*' to 
the adaptive firms model with time lag = 1. 
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Zipf plots of output 5 active regions (not shown here) display the 
same characteristics. 

When the time lag is 5, the most productive region accounts for 
more than 50 perc. of total production. The figure is 18 perc. for 
the second peak. The distribution typically is "winner takes all". The 
equivalent figures when the time lag is 1 are 10 and 8.5 perc. 

In conclusion, the patterns clearly display some intermediate scale 
organisation in active and less active zones: strongly correlated active 
regions are responsible for most part of the production. The relative 
importance of these regions obeys a Zipf distribution. 



3.5 Wealth and production histograms 

The multiplicative random dynamics of capital and the direct obser- 
vation of wealth and production would lead us to predict a scale free 
distribution 6 of wealth and production. 

The cumulative distribution functions (cdf ) of wealth and produc- 
tion observed on figure 8 are indeed wide range and do not display 
any characteristic scale: The data wealth and production were taken 
for the same conditions as the previous figures at the end of the simu- 
lation, i.e. after 5000 time steps. The medium range of the cdf when 
time lag is 1 (figure 8a) extends on one and a half decade with an 
apparent slope of 1 ± 0.05 in log-log scale. 

This observed dependence of the wealth cdf, log normal at lower 
A values followed by power law at intermediate A values, is consistent 
with expressions derived for pdf in the literature on coupled differential 
equations with multiplicative noise. Bouchaud and Mezard 14 e.g. 
obtained: 

P( w ) = Z eXP ~~^ (8) 

(where w stands for the wealth relative to average wealth A), from 
the differential system: 

fJA . 

-j± = m {t)-A i + J-{A-A i ). (9) 



D rather than vertically integrating production, we applyed the clipping, horizontal in- 
tegration and ordering algorithm to firms at the output layer (k = 0) 

6 What we mean here by scale free is that no characteristic scale is readily apparent from 
the distribution as opposed for instance to gaussian distributions. Power law distributions 
are scale free. A first discussion of power law distributions generated by multiplicative 
processes appeared in [TT5] . 
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where rji(t) is a random multiplicative noise, with variance a 2 ; \i = 
1 + 4- 

At higher wealth, the straight line giggles and drops much faster: 
this is because of the underlying region structure. The last 80 perc. 
of the wealth is concentrated in two rich regions and its distribution 
is dominated by local diffusion phenomena in these regions. 

The departure form the standard (equ.8) distribution is even more 
noticeable when avalanches are present. The large wealth shoulder 
is bigger (95 perc. of production) and the first point at zero wealth 
stands well above the rest of the distribution: it corresponds to those 
50 perc. of the firms which are momentarily bankrupted. The fraction 
of bankrupted firms fluctuates in time and so does the slope of the 
linear segment 7 . 




Figure 9: Cumulative distribution of wealth (red '+') after 5000 iteration 
steps. Parameter choices are the same as the previous figures. 

In conclusion, the observed statistics largely reflect the underlying 
region structure: at intermediate levels of wealth, the different wealth 
peaks overlap (in wealth, not in space!): we then observed a smooth 
cdf obeying equation 8. At the large wealth extreme the fine structure 
of peaks is revealed. 

7 both fluctuations are correlated since the slope of the linear segment depends upon 
the number of firms in the distribution 
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4 Conclusions 



The simple model of production networks that we proposed presents 
some remarkable properties: 

• Scale free distributions of wealth and production. 

• Large spatial distribution of wealth and production. 

• A few active regions are responsible for most production. 

• Avalanches of bankruptcies occur for larger values of the time 
lag between bankruptcy and firm re-birth. But even when most 
firms are bankrupted, the global economy is little perturbed. 

Are these properties generic to a large class of models? we will first 
briefly report on equations which display similar behaviour and then 
examine the results which we obtained with variants of the model. 

4.1 Formal approaches of similar dynamics 

A number of models which display equivalent phenomena have been 
proposed and formally solved. We kept our own notation to display 
similarities: 

• Growth by deposition on surfaces |15j. Edwards/ Wilkinson: 

dA 

— = D -AA + rj(x,t) (10) 

A stands for the distance to the interface. D is the surface diffu- 
sion constant of the deposited material and rji(t) is an addititive 
noise. Other models were proposed by Karkar/Parisi/Zhang, 
Derrida / Spohn [T5] . etc. 

• Generalised Volterra-Lotka from econophysics: (Bouchaud|14j. 
Cont, Mezard, Sornette, Solomon 16 etc.) 

rlA ■ 

( ^± = A i . m {t) + Y. J H A o-ll J 3i A i (11) 

j j 

A stands for individual wealth of agents and rji(t) is a multi- 
plicative noise. Agents are involved in binary transactions of 
"intensity" J^j. Mean field formal solutions displays scale free 
distribution of wealth. Simulations display patterns on lattice 
structures (Souma etal|17j). 
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• Solomon etal 18 . Reaction-Diffusion AB models. 
dA 

— = k-A-r](x,t) + D-AA (12) 

A is the chemical concentration of a product involved in an auto- 
catalytic chemical reaction, D is its diffusion constant. Simula- 
tions and formal derivations yield spatio-temporal patterns sim- 
ilar to ours. 



4.2 Variants of the original model 

We started checking three variants, with for instance more realistic 
production costs taking into account: 

• Influence of capital inertia: production costs don't instantly 
readjust to orders; capital and labour have some inertia which 
we modeled by writing that productions costs are a maximum 
function of actual costs and costs at the previous period. 

• Influence of the cost of credit: production failures increase credit 
rates. 

The preliminary simulations confirm the genericity of our results. 

The third variant is a model with "adaptive firms". The lattice 
connection structure supposes a passive reactive behaviour of firms. 
But if a firm is consistently delivering less than the orders it receives, 
its customers should order less from it and look for alternative sup- 
pliers. Such adaptive behaviour leading to an evolutive connection 
structure would be more realistic. 

We then also checked an adaptive version of the model by writing 
that orders of firm i are proportional to the production capacity A of 
the upstream firms connected to firm i. Simulations gave qualitative 
results similar to those obtained with fixed structures. 

We observe that adaptation strongly re-enforce the local structure 
of the economy. The general picture is the same scale free distribution 
of production and wealth with metastable patterns. Due to the strong 
local character of the economy: 

• Avalanches of production are observed (see figure 9), even when 
time lag is short (time lag of 1). 

• The spatial periodicity of the active zones is increased (see figure 
9 with larger density of smaller zones). But again the activity 
distribution among zones is like "winner takes all" (figure 7). 
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Figure 10: Wealth and production patterns for a network of "adaptive" firms. 
The conventions and parameters are the same as for figures 3, 4 and 5, for a 
200x10 lattice. Time lag is 1, the two upper patterns correspond to t = 1500, 
the lower ones were taken when t = 1998. 



4.3 Checking stylised facts 

Even though the present model is quite primitive 8 it is still tempting 
to draw some conclusions that could apply to real economies. The 
most striking result to our mind is the strong and relatively stable 
spatial disparities that it yields. Let us compare this prediction to 
the economic situation of developing countries: large and persistent 
disparities in wealth and production level as compared to developed 
countries. We can even go further and raise questions about the influ- 
ence of the depth of the production network or the kind of investment 
needed: 

• One generally agrees that disparities between developing and de- 
veloped countries increased since industrial revolution. This is 
also a period during which production became more specialised, 

8 We e.g. discuss a "Mickey Mouse" economy with fixed prices independent from supply 
and demand. Introducing price dynamics is not a major challenge: we would simply face 
an "embarras de richesse" having to choose among either local or global prices. In fact both 
kind of adjustment have already been tested: global adjustment in the case of production 
cost connected to production failure through credit costs, or local adjustment in the case of 
adaptive behaviour. We have already shown that they don't change the generic properties 
of the dynamics. 
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which translates in our model as increasing the network depth: 
for instance a shoemaker would in the past make and sell a pair of 
shoes from leather obtained from a cattle breeder. Nowadays the 
shoe production and delivery process involve many more stages. 
Our simulations have shown that increasing depth increases the 
fragility of economies to failures and bankruptcies. The new in- 
dustrial organisation may have detrimental effects on developing 
economies. 

• Obviously investment policies in developing countries yield some 
coordination across the whole production chain. Bad economic 
results might be due to very local conditions but can also reflect 
the lack of suppliers /producers connections. 

The above remarks are not counter-intuitive and these conclusions 
could have been reached by verbal analysis. What is brought by the 
model is the dramatic and persistent consequences of such apparently 
trivial details. 
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